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Abstract 

A report is made on the iuminosity and puise-period evoiution of the Be binary X-ray puisar, GX 
304-1, during a series of outbursts from 2009 to 2013 observed by MAXi/GSC, RXTE/PCA, 
and Fermi/GBM. in totai, tweive outbursts repeated by ~ 132.2 days were observed, which 
is consistent with the X-ray periodicity of this object observed in the 1970s. These 12 out¬ 
bursts, together with those in the 1970s, were found to aii recur with a weii defined period of 
132.189±0.02 d, which can be identified with the orbitai period. The puise period at -- 275 
s, obtained from the RXTE/PCA and Fermi/GBM data, apparentiy exhibited a periodic mod- 
uiation synchronized with the outburst period, suggesting the puisar orbitai motion, which is 
superposed on a secuiar spin-up trend throughout the entire active phase. The observed 
puise-period changes were successfuiiy represented by a modei composed of the binary or¬ 
bitai moduiation and puisar spin up caused by mass accretion through an accretion disk. The 
orbitai eiements obtained from the best-fit modei, inciuding the projected orbitai semi-major 
axis ax sin/ ~ 500 - 600 iight-s and an eccentricity e ~ 0.5, are typicai of Be binary X-ray puisars. 

Key words: puisars: individuai (GX 304-1) — stars: neutron — X-rays: binaries 


1 Introduction 

X-ray binary pulsars (XBPs) are systems consisting of mag¬ 
netized neutron stars and mass-donating stellar companions. 
According to the type of the companion, they are classified 
into several subgroups, including Super Giant XBPs and Be 
XBPs as major members (e.g. Reig 2011). Since these neu¬ 
tron stars are strongly magnetized, the matter flows from the 
companion are dominated by the magnetic pressure inside the 
Alfven radius, and are then funneled onto the magnetic poles 


along the magnetic field lines. Since the accreting matter trans¬ 
fers its angular momentum at the Alfven radius to the neutron 
star, the spin-up rate and the mass accretion rate (i.e. the X- 
ray luminosity) of an XBP are thought to be closely correlated 
(e.g. Ghosh & Lamh 1979). The issue is relevant to the evo¬ 
lution scenario of XBPs and their magnetic fields, and thus 
very important. Although several observational results to ex¬ 
amine the relation have been reported so far (e.g. Finger et 
al. 1996; Reynolds et al. 1996; Bildsten et al. 1997), the ques- 
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tion is still unsettled because we need to select an object with 
a well-determined magnetic-field strength and large luminosity 
swings. Furthermore, we need to disentangle intrinsic period 
changes from orbital Doppler effects. 

GX 304— 1 is a Be XBP discovered by X-ray balloon obser¬ 
vations in 1967 (e.g. McClintock et al. 1971). The source shows 
X-ray properties typical of Be/X-ray binaries; large flux varia¬ 
tions (Ricker et al. 1973), the 272 ±0.1-s coherent pulsation 
(Huckle et al. 1977; McClintock et al. 1977), a hard X-ray spec¬ 
trum represented by a power-law with an absorption column 
density A^h ~ 1 x 10^^ cm“^ and a photon index F ~2 up to 40 
keV (White et al. 1983), and periodic flaring events, namely out¬ 
bursts, separated by ~ 132.5 days presumably attributed to the 
binary orbital period (Priedhorsky & Terrell 1983). The optical 
counterpart was identified with a Be star (Mason et al. 1978), 
which is characterized by double-peak FIq emission, strong Hei 
and Oi absorptions (Thomas et al. 1979; Parkes et al. 1980), 
and photometric variability (Menzies et al. 1981). Based on 
these line features, Corbet et al. (1986) suggested that the cir- 
cumstellar disk is edge-on (i ~ 90°). From the visual extension 
(Av = 6.9 mag.), the source distance is estimated to be 2.4±0.5 
kpc (Parkes et al. 1980). However, the binary orbital elements 
have not been determined, yet. 

Since the last outburst was observed in the 1980s (Pietsch et 
al. 1986), the source had been in a quiescent state with no sig¬ 
nificant X-ray emission for 28 years. Its long absence was bro¬ 
ken by the hard X-ray detection with INTEGRAL in 2008 lune 
(Manousakis et al. 2008). Since then, continuous X-ray mon¬ 
itoring with the MAXI GSC and the Swift BAT has detected 
nine outbursts every ~ 132.5 d, from 2009 November to 2012 
November (Yamamoto et al. 2009; Krimm et al. 2010; Mihara et 
al. 2010a; Nakajima et al. 2010; Ktihnel et al. 2010; Yamamoto 
et al. 2011a; Yamamoto et al. 2011b; Nakajima et al. 2012). 
Across these multiple outbursts, the 2-20 keV MAXI GSC 
intensity varied nearly by two orders of magnitude. Suzaku, 
RXTE and INTEGRAL observations, performed at a bright out¬ 
burst phase, revealed a cyclotron resonance scattering feature 
(CRSF) at around 54 keV (Yamamoto et al. 2011c; Klochkov 
et al. 2012b). The estimated surface magnetic field, Bb — 
4.7 X lO’^^ (1 -I- zg) G where zg represents the gravitational red- 
shift, is almost the highest among those of XBPs of which the 
CRSF signature is significantly detected. In the 2010 August 
outburst, the barycentric pulsation period was ~275.37-275.46 
s, which suggests that the pulsar was spun down by ~ 3 s dur¬ 
ing the 28 years of quiescence (Devasia et al. 2011; Yamamoto 
et al. 2011c). The timing analysis of RXTE data revealed the 
complex pulse-profile dependence on the energy band as well 
as the luminosity, and the quasi-periodic oscillation at ~ 0.12 
Hz (Devasia et al. 2011). The period and pulsed flux have also 
been monifored by the Fermi GBM pulsar project (Finger & 
lenke 2013). All these data thus provide us with a valuable op¬ 


portunity to study the relation between the luminosity and spin 
period of the XBP. 

In this paper, we analyze the X-ray outburst light curve of 
GX 304— 1 during the active period from 2009 to 2013 obtained 
with the MAXI GSC, and the pulse period variations derived 
from the RXTE/PCA and Fermi/GBM data. Our goal is to sep¬ 
arate the orbital doppler effects and the intrinsic pulse-period 
change. The observation and data reduction are described in 
section 2, while the data analysis and results in section 3. We 
discuss the results in section 4. 

2 Observation 

2.1 MAXI monitoring 

Since the MAXI (Monitor of All-sky X-ray Image; Matsuoka 
et al. 2009) experiment onboard the ISS (International Space 
Station) started its operation in 2009 August, the GSC (Gas Slit 
Camera; Mihara et al. 2011), one of the two MAXI detectors, 
has been scanning almost the whole sky every 92-minute or¬ 
bital cycle in the 2-30 keV band. We utilized archived GSC 
light-curve data for GX 304—1 in 2-A keV, 4-10 keV and 10- 
20 keV bands, which are processed with a standard procedure 
(Sugizaki et al. 2011) by the MAXI team and available from the 
MAXI web site*. Figure 1 top shows the 2-20 keV MAXl/GSC 
light curve from 2009 August to 2013 December in 1-d time bin. 
It clearly reveals the recurrent outburst activity by a ~ 132.2 d 
interval, which is consistent with the 132.5 ± 0.4 d periodicity 
found in the 1971-1972 outbursts with Vela 5B (Priedhorsky & 
Terrell 1983). Figure 1 covers 12 epochs of the outbursts pre¬ 
dicted from the 132.2-d period cycle, which we consecutively 
name A, B, C,..., and L. The outburst 1 in 2012 August reached 
~ 2.4 photons cm^ s“* (~ 0.6 Crab) in the 2-20 keV band, 
which is the highest among the flaring events that have ever 
been observed from this source. 

If the source distance of D = 2.4 kpc, the energy spectrum of 
a cutoff power-law with a photon index F = 0.35 and an e-fold 
energy Sfoia = H keV (Yamamoto et al. 2011c) are employed, 
the observed flux of 1 photon cm“^ s“* in the 2-20 keV band 
corresponds to the bolometric luminosity of 1.38 x 10^^ erg s“* 
in an isotropic emission source. We use the conversion factor to 
estimate the luminosity form the MAXl/GSC data, hereafter. 
The scale of the estimated luminosity is shown at the right-hand 
ordinate of the light-curve (figure 1 top). 

2.2 Fermi GBM data 

The GBM (Gamma-ray Burst Monitor; Meegan et al. 2009) 
onboard the Fermi Gamma-Ray Space Telescope is an all-sky 
instrument, sensitive to X-rays and gamma-rays with energies 

' http://maxi.riken.jp 
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Fig. 1. (Top) MAXI GSC 2-20 keV light curve of GX 304-1 from 2009 
August 15 (MJD 55058) to 2014 January 21 (MJD 56678) in 1-d time bin. 
Arrows indicate the expected outburst epochs of the 132.2-d cycie. Tweive 
epochs invoived in the light curve are named A, B, C, .., and L, as labeied. 
Vertical bars in blue indicate the epochs of the RXTE pointing observations. 
(Bottom) Barycentric pulse period (red cross) and pulsed 12-50 keV flux 
(black solid triangle) measured by the Fermi GBM. 

between 8 keV and 40 MeV. The Fermi GBM pulsar project^ 
provides results of timing analysis of a number of positively de¬ 
tected X-ray pulsars, including their pulse periods and pulsed 
fluxes (Finger et al. 2009; Camero-Arranz et al. 2010). We uti¬ 
lized the archived results of GX 304— 1 to systematically inves¬ 
tigate the long-term variations of the pulse timing from 2010 to 
2013. 

Figure 1 (bottom) shows time variations of the barycentric 
pulse period and the pulsed photon flux in the 25-50 keV band 
of GX 304—1, as made publicly available by the project. The 
pulsed emission was detected over bright phases of the nine 
consecutive outbursts from B to J. The period range of 275.2- 
275.6 s during the outburst C is consistent with the results of 
275.37 s from the RXTE data (Devasia et al. 2011) and 275.46 
s from the Suzaku data (Yamamoto et al. 2011c), both obtained 
in the same outbursts. 

2.3 RXTE observations and data reduction 

The RXTE observations of GX 304— 1 were performed on the 
outbursts B, C, D, and E, using the Proportional Counter Array 
(PGA: Jahoda et al. 2006) operating in 3-60 keV and the High- 

^ http://gammaray.nsstc.nasa.gov/gbm/science/pulsars/ 


Energy X-ray Timing Experiment (FlEXTE: Rothschild et al. 
1998) providing 20-250 keV data. As indicated in figure 1 
(top), total 71 observations covered the outbursts C, D, and E 
with a high frequency of almost once per day, each with an ex¬ 
posure of 0.5-18 ks. We analyzed the PGA data to study the 
pulse timing properties with photon statistics and time resolu¬ 
tion better than those available with other instruments. 

Data reduction and analysis were performed using the stan¬ 
dard RXTE analysis tools released as a part of HEASOFT 
6.14, and the CALDB files of version 20111205 provided by 
NASA/GSFC. All the data were first screened with the standard 
selection criteria, that the spacecraft pointing offset should be 
smaller than 0.02, the earth-limb elevation angle be larger than 
10°, and the time since the last SAA passage be longer than 
30 minutes. We used data only from the top layer of the PCU- 

2 unit, which is the best calibrated among all PCUs. In the 
timing analysis, we used the reduced data with Good-Xenon 
data mode. Generic event-data mode (E_125us_64M_0_ls), or 
Generic binned-data mode (B^50ms_128M_0_254), which has 
a time resolution better than 0.25 s. 

3 Analysis 

3.1 Outburst intervals and orbital period 

To quantify the outburst periodicity, we determined their peak 
epochs by fitting each profile with a Gaussian function. Since 
some of the peaks have asymmetric profiles and hence cannot 
be approximated by a single Gaussian function, the fit was per¬ 
formed within a narrow span of ~ 50 d around each peak so that 
the fit is reasonably accepted. If outburst has multiple peaks, 
their peak epochs were determined separately. Table 1 summa¬ 
rizes the obtained epochs, fluxes, and the luminosities estimated 
from the typical energy spectrum at the peaks. 

X-ray outbursts of Be X-ray binaries are known to be largely 
classified into normal-type (type-I) and giant-type (type-II) ones 
(e.g. review by Reig 2011). The normal outbursts emerge near 
the pulsar periastron passage and their peak luminosity reaches 
~ lO®’^ erg s“^, while the giant ones may appear at any orbital 
phase and can be more significantly luminous than the other 
type of events. The regular periodicity of the observed outbursts 
suggests that they are mostly categorized into the normal type. 
Actually, the first seven outbursts from A to G, with a single¬ 
peak profile, satisfy all the normal-type conditions. Similarly, 
the outbursts observed in 1971-1972 with Vela 5B also exhib¬ 
ited the same periodicity and their peak fluxes were in the same 
level (Priedhorsky & Terrell 1983). Therefore, they should be 
categorized into the normal type, too. 

To refine the orbital period, we extrapolated the periodicity 
of the 2009-2012 outbursts (A to G) to that in 1971-1972 as¬ 
suming that they are at the same orbital phase. The top panel 
of figure 2 shows a plot of the peak epochs Tpeak in 2009-2013 
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# of period cycle from MJD 41675.6 (1972 Dec 14) 

Fig. 2. (Top) Epochs of outburst peaks against the number of the 132.2-d 
period cycles counted from that on MJD 41675 (1972 December 14). The 
solid line represents the best-fit linear function to the data of cycle= 0 at MJD 
41675 and cycle=102-108, (Bottom) Residuals of the data from the best-fit 
linear function. The area of the circles are proportional to the outburst-peak 
intensity in both panels. 

against the number n of 132.2-d cycles from the initial epoch 
of MJD 41675.6 (1972 December 14) in Vela 5B era. Thus, 
the epochs Tpeak of n = 0 (MJD 41675.6) and n =102-108 
(outbursts A to G) line up linearly, around which the data of 
n =102-108 scatter by 2.1 d in the standard deviation. We em¬ 
ployed this deviation as the uncertainty in each data, and fitted 
the epochs with a linear function to obtained the best period es¬ 
timate of 132.1885 ± 0.022 and the initial peak epoch of MJD 
41675.0 ± 2.2 with l-rr error. The best-fit function, 

rpeak = 41675.0 +132.189 xn (MJD) (1) 

is drawn together on the data at the top panel of figure 2, and the 
bottom panel shows the data-to-model residuals. Their values 
are listed in table 1. 

As presented in figure 2 and table 1, we also extended the 
best-fit linear function to n = 113 (outburst L). The residuals 
suggest that the outburst periodicity changed after n = 109 (out¬ 
burst H), as investigated below. 

3.2 Evolution of outburst orbital profile 

Figure 3 left shows expanded three-color MAXl/GSC light 
curves of the individual outbursts, as a function of time from 
the peak predicted by equation (1). The hardness ratios are pre¬ 


sented in the right panels. 

The light curves of the outbursts A to G in figure 3 (left) 
reconfirm that their peak epoch agrees, within a scatter of 2.1 d 
with the prediction of the best-fit cycle. However, these profiles 
vary considerably from cycle to cycle. While the outburst D has 
a slight enhancement over the symmetrical profile prior to the 
peak, the outbursts F and G have a subsequent small peak after 
the main peak. 

As expected from figure 1, the outbursts H and I, the latter 
being the brightest, clearly deviate in figure 2 and figure 3 form 
equation (1), because they both show a double-peaked profiles. 
This suggests that these two, or at least outburst 1, are possibly 
the giant outbursts. In the outburst J, no distinct main peak is 
seen and a low-level activity continued for 40 days. 

During bright outburst phases, the two hardness ratios in fig¬ 
ure 3 (right) slightly change with a positive correlation to the 
flux. In addition, the (4-10 keV / 2-A keV) hardness ratio ex¬ 
hibits sharp increases at the peak in the outbursts E, F, and G. 
They correspond to dips in the softest light curve (red), and thus 
can be attributed to an increase of the absorption column den¬ 
sity. 

3.3 Pulse period changes 

To investigate the change in the pulse period with a better pre¬ 
cision and sensitivity than those available with the Fermi GBM 
shown in red in figure 1 (bottom), we analyzed the RXTE PCA 
data, which samples the outbursts C, D, and E on the almost 
daily basis. The data for the outburst C has been already ana¬ 
lyzed by Devasia et al. (2011) and the results revealed a com¬ 
plex pulse profile that depends on the energy and the luminos¬ 
ity. We performed the timing analysis taking account of these 
properties as follows. Figure 4 top and middle panels show the 
MAXl/GSC 2-20 keV light curves in 1-d time bin, and those of 
the 2-6 keV RXTE/PCA count rates per 272 s, approximately 
equal to the pulse period. The large variability in the PCA count 
rate around the outburst peaks represent the variation compo¬ 
nent other than the 275-s pulsation, as Devasia et al. (2011) re¬ 
ported in the outburst C. 

We performed the pulse period search by the folding method 
using efsearch in XRONOS, for the PCA data after the 
barycentric correction was applied with f axbary in FTOOL. To 
optimize the period determination accuracy, the folding analysis 
was carried out by combining data taken within each interval up 
to 2 days. The energy band was limited to 2-6 keV because the 
sharp dip structure in the pulse profile below 8 keV (Devasia et 
al. 2011). The errors on the obtained pulse period was estimated 
by the method of Larsson (1996) using the folded pulse profile 
with 64 phase bins. 

Figure 4 (bottom) show the obtained pulse period against 
the observation epochs, where the results of the Fermi GBM 
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Fig. 3. (Left) Light curves of the outbursts A to L, in 2-4 keV (red), 4-10 keV (biack), and 10-20 keV (green) bands shown as a function of day modulo the 
132.189 d outburst period. Arrows at the top of each panei indicate the peak position determined by Gaussian fit (tabie 1). (Right) Time evoiution of 4-10 keV 
/ 2-A keV (red) and 10-20 keV / 4-10 keV (black) hardness ratios. 
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Table 1. Out burst peaks of the GX 304-1 active period in 2009-2013 observed with MAXI GSC 


OutIDt 

Epoch 

(MJD) 

Residual* 

(d) 

Peak flux^ 

( photons cm^ s“^ ) 

Luminosity’^ 

( 10®’’ erg s"^ ) 

A 

55154.00 ±0.11 

-4.18 

0.198 ±0.005 

0.27 

B 

55289.30 ±0.10 

-1.07 

0.275 ±0.004 

0.38 

C 

55425.52 ±0.02 

2.96 

1.822 ±0.008 

2.51 

D 

55555.70 ±0.03 

0.95 

1.182 ±0.008 

1.63 

E 

55688.66 ±0.03 

1.73 

1.141 ±0.008 

1.57 

F 

55817.98 ±0.03 

-1.14 

1.137 ±0.008 

1.57 

G 

55952.10 ±0.05 

0.79 

1.266 ±0.009 

1.75 

H* 

56075.67 ±0.09 

-7.84 

0.564 ±0.025 

0.78 

H 

56090.93 ±0.21 

7.43 

0.864 ±0.014 

1.19 

r 

56210.72 ±0.05 

-4.97 

0.978 ±0.012 

1.35 

I 

56235.79 ±0.05 

20.10 

1.997 ±0.020 

2.76 

j* 

56337.64 ±0.58 

-10.24 

0.098 ±0.003 

0.13 

j* 

56354.21 ±0.23 

6.33 

0.151 ±0.009 

0.21 

j 

56370.84 ±0.37 

22.96 

0.186 ±0.012 

0.26 

K 

56483.86 ±0.15 

3.80 

0.177 ±0.009 

0.24 

L 

56617.61 ±1.09 

5.35 

0.045 ±0.004 

0.06 


*: Second or third highest peak in a outburst cycle 


Outburst ID designated by the number of period cycles (figure 1) 
Residual of peak epoch from the best-fit period cycle of equation (1) 
Peak photon flux in 2-20 keV band in units of photons cm“^ s“^ 

: Bolometric luminosity estimated from the 2-20 keV photon flux 


pulsar project, the same as in figure 1, are plotted together. 
The RXTE/PCA and the Fermi/GBM results are consistent with 
each other. 

3.4 Determining the binary orbitai eiements 

The pulse-period evolution in figures 1 and 4 is characterized 
by two distinct effects. One is cyclic modulation synchronized 
with the 132.2 d outburst period, considered as due to tbe or¬ 
bital Doppler effects. The other is a secular spin-up trend over 
the 3 years, obviously due to the accretion torque. Thus, the 
observed pulse-period evolution is considered to be a complex 
composite of these two effects. Since both of them are sup¬ 
posed to depend on the orbital phase, it is not straightforward 
to separate them. Hence, we construct a semi-empirical model 
representing the two effects, and then fit it to the data, in an at¬ 
tempt to simultaneously determine the orbital elements and tbe 
luminosity-dependent intrinsic pulse period change. 

Be XBPs are known to exhibit spin-up episodes during 
bright outburst phases (e.g. Bildsten et al. 1997). This is nat¬ 
urally explained by an increase in the accretion rate, and the 
associated transfer of the angular momentum to the neutron star 
via disc-magnetosphare interactions. In the classical theoretical 
model by Ghosh & Lamb (1979), the pulsar spin-up rate — Pgl 
(s yr“^) is given by 

—Pgl = 5.0 x 10 Pso'r^i^a) PspinT/Jy (2) 


where pso, Rnsg, ATns©, hs, Pspin, and L 37 are the magnetic 
dipole moment of the neutron star in units of 10^° G cm^, radius 
in 10® cm, mass in Mq, moment of inertia in 10^® g cm^, spin 
period in s, and the luminosity in 10®^ erg s“^, respectively, 
while njcus) is a dimensionless torque that depends on the fast¬ 
ness parameter ujb, defined by the ratio of the angular velocity of 
the pulsar spin to that of the Kepler rotation of accreting matter 
at around the Alfven radius. 

In the application of equation (2), GX 304— 1 is particularly 
suited for several reasons. First, the value of pso can be ac¬ 
curately estimated from the CRSF detection at 54 keV, which 
implies the surface magnetic field, Bb = 4.7 x 10^^(1 -I- zq) 
G (Yamamoto et al. 2011c). If the canonical neutron-star mass 
THns = I. 4 M 0 and radius Pns = 10® cm are assumed, the grav¬ 
itational redshift zq is 0.3, and then pso is calculated as 

P30~^BbR%s=3.1 (10®°Gcm®) (3) 

(Wasserman & Shapiro 1983). Second, a continuous record 
of L37 is available from tbe MAXI GSC light curve; we cal¬ 
culate it with the conversion factor as described in section 2.1. 
Furtbermore, njws) — 1-4 can be regarded approximately as a 
constant from the long spin period, Pspin — 275 s, such that 
cus -C 1 (Ghosh & Lamb 1979). We also employ Ps calcu¬ 
lated from its approximate equation with Mns and Pns given 
by Ravenball & Pethick (1994) as 
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Fig. 4. Details of the outbursts C (left), D (center), and E (right), (Top) MAXI/GSC light curves in the 2-20 keV band, in 1-d time bin. (Middle) THe 2-6 keV 
RXTE/PCA count rate per 272-s time bin during the pointing observations. (Bottom) Barycentric pulse periods obtained from the RXTE/PCA data (blue dot) 
and the Fermi/GBM pulsar archive data (red cross) All vertical error bars represent the 1-cr uncertainty. 


/45~0.21- 


Mns^^ns 


= 1.0 gcm^). 


1 — 2GMns/J?nsc^ 

Combining these pieces of information, equation (2) yields 


(4) 


-Pgl if 1.7 X (sd-i). 


(5) 


Equation (2) implies that the spin-up rate —Paptn follows 
the luminosity L as — Psptn oc Lr. The relation has been cali¬ 
brated with actual X-ray data (Reynolds et al. 1996; Bildsten et 
al. 1997). Besides this, the comparison of absolute spin-up rates 
with equation (2) has been hampered by a large uncertainty in 
the bolometric luminosity correction, which is in turn due to 
beaming effects (e.g. Bildsten et al. 1997; Wilson et al. 2002). 
We hence employ a spin-up model expressed by —Pspin = 
in which the power-law index and the constant coefficient a are 
treated as free parameters. 

XBPs are also known to spin down during the quiescence 
due to the propeller effects. In fact, GX 304—1 exhibited a spin 
down by ~ 3 s for the 28 years of quiescence from 1980 to 2008, 
with an average period derivative ofP~3x 10“® ss“^. The 
rate, though much smaller than the spin-up rate during bright 
phases ( P ^ —10“^ s s“^), may not be negligible. We ac¬ 
counted its contribution with a constant spin-down parameter, 
/3, added to Ppin as an offset. 


By using the spin-up and spin-down models discussed 
above, the intrinsic pulsar-spin period Pspin(t) is expressed as 

Pspin(f) = Po-f f Pspin(r)dr (6) 

Jtq 

= Po+ / {-aLl,{r)+p}dr (7) 

Jtq 

where we set the time origin tq at the periastron passage in the 
outburst C and denote the pulse period at to as Po = Pspin(To). 

The period modulation due to the binary motion, can be cal¬ 
culated using the binary orbital elements, which consists of the 
orbital period Pb, the eccentricity e, the projected semi-major 
axis flxsini, the argument of the periastron cjo, and the epoch of 
periastron passage which we identified with tq defined above. 
The pulsar orbital velocity vi (t) along the line of sight is given 
as 

, , 27raxSini r / / \ \ i /on 

vi{t) = - |cos(t/(f)-Ea;o)-|-ecosajo| (8) 

PbV 1 - 

where i/{t) is a parameter called “true anomaly” describing the 
motion on the elliptical orbit, and calculated from the Kepler’s 
equation. The observed barycentric period, Pobs(f), is then ex¬ 
pressed by 
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Fig. 5. (i) Bolometric luminosity estimated from the 2-20 MAXi/GSC iight curve in 1-d time bin. (ii) Barycentric puise period measured with the RXTE/PCA 
(biue) and the Fermi/GBM (red), fitted with equations (7)-(9). Biack and green soiid lines represent the best-fit model with free 7 and that with 7 = 6/7, 
respectiveiy. A dashed line represents the period moduiation due to the orbital Doppler effect, (iii) Residuals of the data from the best-fit model with free 7. (iv) 
Residuals from the 7 = 6/7 model. 


Pob.{t) = Pspin(t) (^1 + (^1 - 

=^Pspin(i) (9) 

We applied the model Pobs(i) consisting of equations (7)- 
(9) to the observed barycentric periods obtained from the 
RXTE/PCA timing analysis in section 3.3 and Fermi/GBM 
archival products. The model involves nine free parameters; 
Po, OL, (3, and 7 in equation (7), and orbital elements in equation 
( 8 ). We utilized the luminosity calculated from the MAXI/GSC 
2-20 keV light curve in 1-d time bin, as described in section 
2.1. The errors associated with the L 37 determinations propa¬ 
gate into those of equation (7), but this gives 1-a uncertainties 
only by ~ 0.003 s, which is smaller than those in the individ¬ 


ual pulse-period measurements, ~ 0.01 s. We fixed the orbital 
period Pb at 132.189 d derived from the outburst cycle over 
28 years in section 3.1, and examine the fit with 7 fixted at the 
theoretical value of 6/7 or free. 

The model has been found to represent the data well. Figure 
5 shows the results of those model fits (with 7 = 6/7 and 7 
free), and table 2 summarizes the obtained best-fit parameters. 
The obtained epoch of the periastron of the outburst C, tq ~ 
55425 (MJD), agrees with that at the outburst flux peak (table 
1) within the uncertainty of the peak-epoch periodicity, ~ 2.1 
d. Comparing the two fit results, the fit with 7 free has a 
lower by 108.3 for a decrease of only by one in the degree of 
freedom (dof). In figure 6, pulsar binary orbits suggested in the 
two fitting models, are illustrated. 
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Table 2. Best-fit orbital parameters 


Parameter 

Model 

7 = 6/7 

7 : free 

Hb (d) 

132.189 (fix) 

axsini (It-s ) 

498 ±6 

601 ± 38 

e 

0.524 ±0.007 

0.462 ±0.019 

ro (MID) 

55425.0200 ±0.0010 

55425.6 ±0.5 

Ho ( s ) 

275.4441 ±0.0008 

275.459 ±0.008 

too ( degree ) 

122.5 ±0.4 

130.0 ±4.4 

a/0.017 

0.2989 ±0.0004 

0.2467 + 0.0010 

,3( 10 -®s s-l ) 

3.27±0.02 

< 0.23 

7 

6/7 (fix) 

1.243 + 0.031 

Jm ( Mq ) 

7.6±0.3 

13.3 + 2.5 

X® / dof 

276.3 / 103 

168.0/102 


All errors represent l-a confidence limit of the statistical uncertainty. 



-500 ' ^ ^ 0 ^ ^ ’ 500 


Fig. 6. Pulsar binary orbits obtained in the two period-change models with 
7 free (biack) and 7 = 6/7 (blue), projected on a plane including the line of 
sight. 

4 Discussion 

4.1 Lumonosity - spin-up relation 

As a result of the model fit to the pulse-period evolution, we 
successfully calibrate the relation between the the luminosity L 
and the spin-up rate, —Papin, which is considered to reflect the 
accretion manner. According to the model of Ghosh & Lamb 
(1979), the power-law index 7 becomes 6/7 in a disk accre¬ 
tion and 1 in a wind accretion. The best-fit 7 obtained from 
the model of 7 free was 1.24 ± 0.03, which is significantly 
larger than 6/7 ~ 0.86 in the disk accretion. Meanwhile, the 
present result agrees with those of other typical Be XBPs, A 
0535-f 26 (Finger et al. 1996; Camero-Arranz et al. 2012), EXO 
2030+375 (Reynolds et al. 1996; Wilson et al. 2002), and 2S 
1417-624 (Finger et al. 1996; inam et al. 2004). As Wilson et 
al. (2002) suggested, the result of 7 > 1 may imply the effect 
of either the spin down or the wind accretion during the low 


luminosity periods. In fact, the fit with 7 = 6/7 yielded the 
significant spin-down of /3 = 3.27 ±0.02 x 10~® s s“^ during 
the outburst intermissions. The rate is consistent with the aver¬ 
age period derivative of ~ 3 x 10“® s s“^ during the 1980-2008 
quiescence. 

Comparing the two model-fit results (7 = 6/7 and 7 free), 
goodness of fits (x^ / dof) certainly prefers the one with 7 free. 
However, neither is acceptable statistically. As seen in figure 
5, the difference between the two models is mostly in the out¬ 
burst intermission in which the period data are not available. 
Therefore, it is hard to conclude which model better represents 
the real situation, from the present result. The problem will be 
solved if the pulse period throughout the entire orbital phase is 
obtained. 

From the best-fit parameters, we also obtained the normal¬ 
ization factor of the absolute spln-up rate, a, compared to that 
predicted in the model of Ghosh & Lamb (1979), which is 
0.017 s d“^ as derived in equation (5) in this source. The ra¬ 
tio, a/0.017 ~ 0.3 in both the two model (7 free and 7 = 6/7), 
Is significantly smaller than unity, like in some other Be XBPs 
(Wilson et al. 2002). This discrepancy is reasonably explained 
by the uncertainty in the luminosity estimate. We calculate the 
luminosity from the observed fluxes sampled by MAXI/GSC 
scan survey in 1-d time bin assuming that the source emission is 
isotropic, as described in section 2.1. Because the observed flux 
is pulsating, the source emission is not completely isotropic. 
The obtained result corresponds to the beaming effect by a fac¬ 
tor ~ 3. 

4.2 Binary properties as Be/X-ray binary pulsars 

From the obtained orbital elements, the mass function /m is 
calculated as 

j. (M*sini)® 

(M*±Mns)2 ^ G Pi ’ ^ 

where G is the gravitational constant, and M* Is the mass of the 
companion star. Table 2 yields /m = 7.6 ± 0.3Mq if 7 = 6/7, 
and 13.3 ± 2.5Mq if 7 is left free. Assuming the mass of the 
pulsar to be Mns = IAMq, the condition of sini < 1.0 con¬ 
strains the companion mass as M* > 9.9Mq for /m = 7.6Mq 
and > 15.7Mq for 13.3Mq. This agrees with the typical Be- 
star mass. The system is actually suggested to be edge-on 
(i ~ 90°) from the emission and absorption lines in the optical 
spectra (Corbet et al. 1986). 

The obtained orbital eccentricity of e ~ 0.5 is typical of a 
Be/XBP subgroup that is characterized by a long spin period 
Hapin ^ 100 s and a long orbital period Pb ^ 100 d (Knigge et 
al. 2011). The value of e ~ 0.5 and 73s = 4.7 x 10^^ (1 ± zg) G 
of GX 304— 1 are in line with the positive correlation found by 
Yamamoto et al. (2014) between these two quantities. 
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4.3 Evolution of outburst light-curve profiles and Be 
disk 

The outburst orbital profile is considered to reflect the spatial 
extent of the circumstellar disk around the Be star on the pulsar 
orbit. During the 3.5 years of active period from 2009 October 
to 2012 November, the profile changed through a series of nor¬ 
mal outbursts (A to G), giant outbursts with a large peak lu¬ 
minosity (H and I), and decay phase (J to L). This evolution 
behavior vary much resembles those of other Be X-ray binaries, 
including A 05354-26 (Camero-Arranz et al. 2012; Nakajima et 
al. 2014), and GRO J1008-57 (Kuhnel et al. 2013; Yamamoto et 
al. 2014). Therefore, the same physical mechanism is consid¬ 
ered to be working in those systems. To explain this, a couple 
of scenarios, including the global one-armed oscillation and the 
precession of the Be disk, have been proposed (Nakajima et al. 
2014 and references therein). The detailed modeling is out of 
the scope of the present paper. 

We thank all the ISS/MAXI team members for their ded¬ 
icated work to enable the science mission operation. We 
also thank the Fermi/GBM pulsar project for providing the 
useful results to the public. This research has made use 
of RXTE data obtained from the High Energy Astrophysics 
Science Archive Research Center (HEASARC) provided by 
NASA/GSFC; The present work is partially supported by the 
Ministry of Education, Culture, Sports, Science and Technology 
(MEXT), Grant-in-Aid No. 24340041 and 25400239. 

References 

Bildsten, L., et al. 1997, ApJS, 113, 367 

Camero-Arranz, A., Finger, M. H., Ikhsanov, N. R., Wilson-Hodge, C. A., 
& Beklen, E. 2010, ApJ, 708, 1500 
Camero-Arranz, A., et al. 2012, ApJ, 754, 20 
Corbet, R. H. D., et al. 1986, MNRAS, 221, 961 

Devasia, J., James, M., Paul, B., & Indulekha, K. 2011, MNRAS, 417, 
348 

Finger, M. H., et al. 2009, arXiv:0912.3847 

Finger, M. H., Wilson, R. B., & Chakrabarty, D. 1996, A&AS, 120, 209 
Finger, M. H., Wilson, R. B., & Harmon, B. A. 1996, ApJ, 459, 288 
Finger, M. H., & Jenke, P. A. 2013, The Astron. Telegram, 4812, 1 
Ghosh, P, & Lamb, F. K. 1979, ApJ, 234, 296 

Inam, S. Q., Baykal, A., Matthew Scott, D., Finger, M., & Swank, J. 2004, 
MNRAS, 349, 173 

Huckle, H. E., Mason, K. O., White, N. E., Sanford, P. W., Maraschi, L., 
Tarenghi, M., Tapia, S., 1977, MNRAS, 180, 21P 
Jahoda, K., Markwardt, C. B., Radeva, Y., Rots, A. H., Stark, M. J., 
Swank, J. H., Strohmayer, T. E., & Zhang, W. 2006, ApJS, 163, 401 
Klochkov, D., et al. 2012, A&A, 542, L28 

Knigge, C., Coe, M. J., & Podsiadlowski, P. 2011, Nature, 479, 372 
Krimm, H. A., et al. 2010, The Astron. Telegram, 2538, 1 
Kuhnel, M., et al. 2010, The Astron. Telegram, 3087, 1 
Kuhnel, M., Muller, S., Kreykenbohm, I., et al. 2013, A&A, 555, A95 
Larsson, S. 1996, A&AS, 117, 197 


Manousakis et al., 2008, Astron. Telegram, 1613 

Mason, K. O., Murdin, P. G., Parkes, G. E. Visvanathan, N., 1978, 
MNRAS, 184, 45P 

Matsuoka, M., et al. 2009, PASJ, 61, 999 
Meegan, C., et al. 2009, ApJ, 702, 791 

McClintock, J. E., Ricker, G. R., Lewin, W. H. G., 1971, ApJL, 166, L73 
McClintock, J. E., Rappaport, S. A. Nugent, J. J., Li, E. K., 1977, ApJL, 
216, L15 

Menzies, J., 1981, MNRAS, 195, 67P 

Mihara, T., Makishima, K., & Nagase, F. 2004, ApJ, 610, 390 

Mihara, T., et al. 2010a, The Astron. Telegram., 2779, 1 

Mihara, T., et al. 2010b, The Astron. Telegram., 2796, 1 

Mihara, T., et al. 2011, PASJ, 63, 623 

Nakajima, M., et al. 2010, The Astron. Telegram, 3075, 1 

Nakajima, M., et al. 2012, The Astron. Telegram, 4420, 1 

Nakajima, M., et al. 2014, PASJ, 66, 9 

Parkes, G.E., Murdin, P.G., Mason, K.O., 1980, MNRAS, 190, 537 
Pietsch, E., Collmar, W., Gottwald, M., Kahabka, P, Ogelman, H., 1986, 
A&A, 163, 93 

Priedhorsky, W. C., Terrell, J., 1983, ApJ, 273, 709 
Reig, P.2011, Ap&SS, 332, 1 
Ravenhall, D. G., & Pethick, C. J. 1994, ApJ, 424, 846 
Reynolds, A. P, et al. 1996, A&A, 312, 872 

Ricker, G. R., McClintock, J. E., Gerassimenko, M., & Lewin, W. H. G. 
1973, ApJ, 184, 237 

Rothschild, R. E., et al. 1998, ApJ, 496, 538 
Sugizaki, M. et al. 2011, PASJ, 63, 635 

Thomas, R.M., Morton, D.C., Murdin, P.G., 1979, MNRAS, 188, 19 
Wasserman, I., & Shapiro, S. L. 1983, ApJ, 265, 1036 
White, N., Swank, J., & Holt, S. S., 1983, ApJ, 270, 711 
Wilson, C. A., Finger, M. H., Coe, M. J., Laycock, S., & Fabregat, J. 
2002, ApJ, 570, 287 

Yamamoto, T., et al. 2009, The Astron. Telegram, 2297, 1 
Yamamoto, T., et al. 2011, The Astron. Telegram, 3309, 1 
Yamamoto, T., et al. 2011, The Astron. Telegram, 3624, 1 
Yamamoto, T., et al. 2011, PASJ, 63, 751 
Yamamoto, T, et al. 2012, The Astron. Telegram, 3856, 1 
Yamamoto, T, et al. 2014, PASJ, 66, 59 



